f=inline('(x-1)/(x+1025)');
x0=1;
for i=1:20
    x0=f(x0);
    fprintf('i= %i,x0= %f\n',i,x0);
end